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MRF.OIL PROFILE IN A NON-UNIFORM FLOW 


Jan Polasek* 


/ 533 ** 

We will discuss a method for calculating the flow around a profile 
in a' nonhomogeneous flow field. This method is based on the idea that 
the profile is replaced by a continuous distribution of vortices and 
sources along the profile skeleton,. • The method is suitable for solving 
the direct and the indirect problem - of profile theory. 


The results given can .be applied for" calculating the air foil in an 
homogeneous and a non-homogeneous flow field, or for calculating 
straight and radial turbine blades, and guide blades, etc. For prac- 
tical applications, the results are: given in the form of formulas, 
suitable for numerical calculations. Examples are also given. 


SYMBOLS : 
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y 8 

t 

y mt 
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x. 


P J y p 


R : 

R. 




profile skeleton chord- 
profile ‘skeleton ordinate 

distribution of profile thickness (measured from skeleton) 
ordinate of profile Central line 

symmetric thickness distribution (measured from the cen- 
tral line) 

coordinates of points on the profile 
curvature radius of leading edge 
curvature radius of trailing edge 

trigonometric auxiliary variable (identical with the 
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polar angle in the complex auxiliary plane it,) 
velocity components 
contour velocity 
basic velocity 

total circulation over profile /6l8 

circulation distribution 
source distribution • 

coefficients in the expansion of profile thickness 

coefficient in the expansion of skeleton direction factor 

coefficient of the expansion of the directional factor of 
central line 

coefficient in the expansion of circulation distribution 

coefficient in the expansion of source distribution 

coefficients in the expansion of components of primary 
velocity 


COMPLEX PLANE Z 

The same notation is used in the physical plane, but the quantities 
are considered dimensionless. The lengths are referred to one-half of 
the chord length (c/2) and the velocities are referred to the basic 
velocity (U q ). 

"iu complex velocity in plane z 

Uy, v velocity components induced by vortices 

u , v velocity components induced by sources 

nl M. 

u K , v K components of contour velocity 


COMPLEX PLANE Z: 


X, Y 
X v> Y v 

m 

U, V 
U Y , V Y 

u , v 

q 5 q 
d k> v k 


rectangular coordinates Z = X + iY 
coordinates of points on the profile 
complex velocity 

components of primary velocity in plane Z 
velocity components induced by vortices 
velocity components induced by sources 
components of contour velocity 


k.'\» 
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COMPLEX PLANE g: 


P, ^ 
£ 

V 


Polar coordinates :% ~ 

distance of point from unit circle (e =p - 1) 

complex velocity in ? plane 

radial and azimuthal velocity components 


INTRODUCTION 

The idea of Birnbaum of using a’’ continuous vortex [distri bution to 
replace the flow around a thin profile has been found to be- ver'yufruit- 
■ ful. The original paper of Birnbaum [T] was followed by j similar j . o. 
[studies? [2-6] , and several applications were considered: Prom the flow 

around a thin profile in a homogeneous flow to the flow around a pro- 
file in a nonhomogeneous or non-steady flow field. Then this was ap- 
plied to the solution of guide blades, 'and blade cascades. Another ad- 
vance was the consideration of the'blade thickness, which can be 
. achieved by a suitable distribution of sources .on the profile skeleton. 

/6l9 

In the calculation of the flow around profiles with a moderate cur- 

i 

vature, a further simplification is introduced by decomposing the vor- 
tices and sources along the prof ile . chord, instead of along its skele- 
;ton. The -calculations become cleared with these simplific ations, and 
they can be performed -'without great mathematical knowledge . The so-called 
‘first-order theory gives very good results for the ratio between the 

shape of an infinite thin profile and the circulation distribution. How- 

) ' 1 

ever, the calculation of the contour velocity and the influence of the 

finite thickness, especially near the leading edge., lead, to certain 

difficulties and inaccuracies. Various authors have tried to eliminate 

1 

or reduce this by introducing various correction factors. Deeper mathe- 
matical investigations have shown that in the derivation of the rela- 
tionships between the shape of an infinite thin profile and the circu- 
lation distribution in a homogeneous flow field, only the "third order" 
terms were ignored, because the "second order" terms drop out automati- 
cally (as can be seen by comparison with the results of paper [6], for 
the case of a moderately curved profile in a homogeneous flow field) . 
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In other eases, especially when calculating the velocity distribution 
along the blade contour and the influence of finite profile thickness, 
the ignored "-second order" terms become noticeable and have a substan- 
tial influence, especially near the leading edge. 

For these reasons, in this paper we developed a complete "second 
order" theory for- the flow around a profile which has moderate cur- 
vature and is not too thick ; in a nonhomogeneous flow, field.. The 

paper is given in two parts: in the first part, we derive the theory 

* 

and in the second part we -give the theoretical results in terms of 
formulas for numerical calculations. These formulas are then used for- 
the practical solution of several selected- problems. First, we will 
calculate the flow around a symmetric and curved profile in a homo- 
geneous flow field, and then we will consider the flow around the same 
profile in a grid arrangement. 


- FIRST PART 


ORIGINAL PAGE IS 
OF POOR QUALITY! 


THEORETICAL ANALYSIS 

2. Infinitesimally-thin profile in a- non-homogeneous flow field 

i 

The profile under consideration is placed in the complex plane z = 
x + iy in such a way that one chord lies along the x-axis, and the ori- 
gin of the coordinate system is at the center. The lengths are con- 
sidered dimensionless and are referred to one-half of the chord length 
(c/2) . 


The profile equation is: 


?/= sS-'cp SpL? 


P 


1 4 Vv. i 


■U'.-i 


where 




( 2 . 1 ) 


( 2 . 2 ) 


/620 


T 


Instead of the variable x, we will; introduce the trigonometric auxi- 
liary variable & with the equation ( • 
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( 2 . 3 ) 


# =2 _ COS & j 0 ^ a- <g 7 T 

For the calculations, it seems appropriate to express the differen- 
tial quotient of the function y (x) in the form of a trigonometric 

s 

series : 


du 00 

-s *.*»**• 


(2.4) 


n=0 


Figure 1 



The profile shape is then given by the expression 

y. ■= ” cos,<> ) + 

CO 

■ + E- 




‘ ^41 (1 — cos nfl-J 

n << 


(2.5) 


The coefficients must satisfy the following condition because of 


equation (2.2): 


S 


b 2m 


i M 2 — 1 . 


( 2 . 6 ) 


Also the components of the velocity of the primary inhomogeneous flow 
field at the profile location, as well as the circulation density (re- 
ferred to the length unit of the chord), are expressed as a continuous 
distribution of vortex threads on the profile using trigonometric series 


u = 1 ,-b Yi it. cos 

<xy 

V ~ £ C0S 
n=o 


(2.7) 


( 2 . 8 ) 
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(2.9) 



The velocitie s and circulation density are considered dimension- 
less, and are referred to a suitable basic velocity U . The coeffi- 
cients in the expansions (2. -4), (2.7), (2.8), and (2.9) are assumed to 
be small of first order, in the sense that the quantities p n , v n , g n 
and are of first order compared with 1. For large n, they go '-to 
zero like the terms of a geometric series with a quotient smaller than 
1. The theory is developed as a second order theory in such a way that 
everywhere the produ cts' of two first-order variables are considered as 
second-order variables, and the products of three first-order variables 
are ignored. 

In practice, it is sufficient to approximate the trigonometric expan- 
sions (2.4) , (2.7), and (2.8) using trigonometric polynomials having a 
relatively low order. 



\y 

■ r ' 


z-x-iy. 

Z~X+iY 

-/ 1 

0 - 1 x -i 

0 / X 

Figure 2 

' 1 

• 


In addition, we will introduce the complex plane Z = X + iY and the 
complex function: 

3 + */(#) (- 2 . 10 ) 

It maps the plane Z onto the plane z in such a way that the distance 
-1 <_ X <_ 1 along the real axis of plane Z is mapped onto a curve y in 

the plane z. In the plane, Z, we introduce for the variable X in the 

interval -1 < X <_ 1, Y=0a trigonometric auxiliary variable & as 

_ ^ , , Y * i 

The circulation density on the profile is ^ ri + (HuJivV f 1 * ' * The Quan- 

t-ities y and dy s /dx are subjected to the restriction that the differ- 
ence between y s and y is small to third order, and therefore, we can 

set y s - y for our theory. 
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follows : 


X = — tos » , o g 


( 2 . 11 ) 


The function f(Z), for example, can be written in the form 

/(£)=< ' (2-12) 

, * > * 1 

The following relationships can hold between the coefficients g n and /^22 
the coefficients B n in the expansion of the directional factor (2.4), 
which are linear relationships 

■ (4 = 4, - 4 + fc| 

• * (ft» - P.J ~ 4 - 34, + J5i? 5 - (2.13) 

-r d (Pj -p 3 )-lJS a - 16174 +-.w, < 

: ; 1,, i « . '• « i 1 • i i > ; ) . ) . i . 

t . ,J ►, V ...... . 

Since we will perform our analysis in the image plane Z, it is ne- 
cessary to map the velocity of the primary flow field from the plane 
z onto plane Z. For the mapping of the primary velocity, we have 


-fr = U — it = (4 



(2.14) 


We have. the following relationship for the derivative of the mapping 
function along -1 <_ X <_ 1, Y = 0: 


dz __ . 


+ i 


2 



(2.15) 


so that the expressions for the components of mapping of the primary 
velocity* along the interval -1' <_ -X <_ 1 along the real axis of the plane 
Z are given by the following: 


# = n -j- <o v ' 

’ <1* (2.16) 

d& ; 


*The mapping (2.10) maps the flow field in the plane z again onto any 
flow field in plane Z - which is a mapping of the flow field from plane 
z. Expressed more simply, this mapping flow field is simply called a 
flow field in the plane z. 
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After substituting from (2.4), (2.?), and (2.8), we find: 

; -fr == A + f*. m. cos \ X - f v„ b 6 -i- J; vi + X § [ y, tii-k -h; 

f. ,, i‘ _ «■’».*; A ) t—o J ’ 2'hTxlto. . J 

< '* OP _ * 1 T 


+ r SK 

i'-O . 


cos-id)-, 


K*o». 44 - “ L + ,g, 14. X - 4£X ^ 4 4 

; s ■ -»_’0 ■ ^ V , .no / i ii-t; n d 


(2.17) 


tf 


' .• .■■J r ((4 $»'+*. f, &!.■) 

t - J 


COS’. 


/ 

The complex velocity in the plane z, which is induced at the profile 
(2.5) by the distributed vortices (2.9), is .given by: 


7623 


it 


~ ^ -T " WZ\ . " 7 " r" > 


( 2 . 18 T 


With the mapping (.2.12), the velocity ' (2 .18) is transformed to 

; b r - i v y = fi 4- ;M \— (' — j ( 219 ) 

: 1 ' az\ 2tc: j- j X - Xjj -P j [Y 4 - f -/ Ujj] ! * 


Along the segment -1 < X < 1 of the X axis, we have: 

i — V'nVHTV’ . 

•- " ! d .(»: 7 2 1 A — A " 0 -f «(//„,— t( a o). « 

The expression (2.20) is transformed into 


, 7 j 

t y. •' • *1® : dHitL,:?- <%; . „• u-iu i« •<$•* v 


VI* + 

: . :V . 


X-A’ 0 ".iltc'l 


and the integrand is expanded according to powers of 




(2.20) 


(2.21) 


< tda;] - X b . : <U if j 

t ’ •' 1 1 V V ; i 

If we restrict ourselves to second-order terms in this expression, 
then we can write 


' t r - ? Vy — — ( l fl + .~ V -'-. y A 

2u.'_,L dili* X-Xo ). 


%(Xd «X i I 

* 14 

• ■ ;t - x 0 . ] 


or, after substituting the variable??,, according to (2.11), we have: 
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( 2 . 22 ) 


U y - i V. ( = ±VU + if JXcos ?i0 — -i- f 7>'„ + 
, ; 0 [ L*-0 *- V 

+ y -fin - i — eOS ti-X " COB-W fr 

ntj -»■ COS X “ C °S ^ 


2<?o (1'4-cos X) -f .</„ [cos (u — i) — X -cos (« -h 1) x] } — - 

J cosy. 


dx 


COS5> ' 


The incomplete integrals on the right side of (2.22) are to be con- 
sidered in terms of the Cauchy principal value. For example, they were 
calculated in [6], Appendices I and IV. Therefore, we will only give /624 
the final result: 


V y — iVy 


~ ± jffo C 0 tg ~ -f 5] g n sinji&j + i |# 0 — g n boS^-i)j -f 

If; y , , - 

Oa "l" ffijj 4- — -#2 |{/o + ffa) + — B a Hr -j- ifr j -{- ■ 


+ 


7 i , 

+ — Ho + . . . -f 


H-tosS. 


({«.,+ f ».] 
+ b. 

4- cos 29-. I 

>1 


3-oos 3 V. 

• 

{ 8 i \ 




i i 

9b + jg-ftj + %5 9o + • • 


4 * 


+ 


£ 

3 




5 

~3 


(2.23) 


The + sign applies for the top side of the interval -1 < X < 1, and 
the - sign for the lower side of this interval. In order not to have 
to write both signs for the velocity, we will consider the angle in 
the interval [-ir, tt ]. The positive value of angle & will refer to the 
top side, and the negative value will refer to the lower side of the 
interval (Figure 3). 


Figure 3 


J. * .&> o 


-/ 




Y 




X 

ti-K < 
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Since the segment -1 <_ X 1 lies along the stream line, then the 
Y-th velocity component must be 0 there, and therefore 

7.4- Vy — 0. . (2.24) 

If we substitute expressions (2.17) and (2.23) for V and V in equa- 
tion (2.24), and if we compare coefficients for equal multiples of the 
cosine of the angle ^ then we find a system of equations which gives /625 
the relationship between the profile shape and the circulation distri- 
bution : 

= — #o + v o — s*o B o — — B hi 

i « «, ' (2.25) 

<J n = B n — -f — Yi (** B n-k + — £ B k)j n— 1,2 } . : 

f k~0 -4 

It should be realized that system (2.25) also results from system 
(2.11) in paper [6], if we set w = 0 there. 


The contour velocity on the profile in the z-plane is obtained from 
the X-th velocity component along the real axis in the plane Z, when 
this velocity component is divided by the modulus of the mapping func- 
tion: 


[i 4- (dy,/a®)»r 


(2.26) 


If we substitute expressions (2.4), (2.17), and (2.13) in equation 
( 2 . 26 ), and if we restrict ourselves only to second-order terms, we 
find: 


10 „ = 1 Yi Pn COS 4- 


9b cotg—8-4- 


eo » * 

£ g n sin tel- -f- 


4- ~fv 0 Bo 4- S V, i? fc ] + ~ fj[ £ v* J B„_* 4- S (vA-U- 4- V, ,i3 fc )1cos «$■ 

/ -4 n~lLt»0 J.=»0 J 

- + t - T s It £ B k B »-h + £ B k + 

* V J “ n-i L “ fc-« \ fe-0 j 

+ [I ' B ‘ ["’ + i Ji ) + i *• (*■ + i *) + t A >{* + t 9, j + . 

+ ~ lh !7° 4- — B h g 0 -f . . . j + 

+ [ B «(|-9. +|- «>)+ Jfc) + + 


(2.27) 


(continued) 
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(continued from previous page) 


4- 


4- — ffi) + B s g 0 + « • * j COS $ + 

[^3 + “0i ) + 4'-“ J? «ft 4 ..i C0S2$4 

| |-g- So 4- ~ <Ji j A- ^ # 6 ffa -4 • • ■ j COS 3$ 4 


+ 


4" f B& ‘Oo 4 • * • j cos 4-D- 4* * • » . 


L3 


. oc 


(2.27) 


/62 6 


It should be realized that the expression (2.24) also follows from 
expression (4.5) in [6], if we set w = 0. 


3. Velocity Field of Sources Distributed on the Skeleton 

The finite profile thickness is considered by decomposing the vor- 
tices and the sources on the profile skeleton. The density of the con- 
tinuous distribution of these sources is considered in the form of 
trigonometric series, just like a circulation distribution: 


« =144 (dy e /da;) 2 ] Vl q, = 2 (q 0 oolg ~ 9- 4 

* , i 2i 

1 °° 

4* q 0 tg— «■ 4 Yi ffi» sin »&), 0 ^ ^ tc. 

2 n-l , ' 


(3-D 


The first term 2 2o'cotg ^ considers the rounding of the leading edge 

1 “ J -’- - ' ' We have 


and the second term 2^ tg — 8,i the rounding of the trailing edge. 

2 

= 0 for the profile with a sharp trailing edge. The total yield of 


‘o 


all of the sources distributed on the skeleton is given by 


jj q, ds = (j rjdcc = 2tt j^. 


ffo 4- 2o 4 — 

2j 


(3-2) 


Figure 4 



original ? age ® 
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Since the profile is closed, the total yield of the sources must be 
zero, which results in a condition which the first three coefficients 
in expansion (3-1) must satisfy: 

& + + ■ (3 - 3) 

* a 


The complex velocity induced by the vortices on the skeleton is, 
found from the velocity induced by the vortices, and this is done by /627 
replacing the coefficients g by r the_ coefficients^^ . The entire ex- 


n 


pression is multiplied by the factor -i, (calculation of the term with 

;$g~& , which does not exist in a circulation distribution, is similar 
2 . f r* 

to the calculation of the term with eolg-r-j)- ) . After mapping into the 

i. . 2_'-i 

Z plane, the velocity on the segment -1 < X < 1, Y = 0 is given by the 


expression: 


- iV, ~ q 0 ~ ffrt ~ S L cos m3 ~i fg 0 t-oCg 1 — S> -j- 
+ 2o 4' 4^ + V q H xiuv‘d\ - 


/ 1 


‘l | |^i + + ila + ” ?ij + 

j \ qAj — — ^3?i + • • 1 


+ — ( + — JS 4 -f . . . j I g ( , ~ q\ t + 

r> 


+ 


l ±Bt + , . .) + 4o + ?J ‘) + j^ a + 2* s ) £ 

(I s - 


+ 


%+-r ( h - 4 + • • - 


COS 3 + 


<h - 

3 + 


+ ~ #5 + • • • j jso + §o 4- — 2ij + 


+ ~ ih (<Zo — 2o) "I" ■ 
o 


cos 2 t) -j~ 


I- f /n + — <h -1 

o j 


I —14 Cg ( i “ 2o) {-••• 

o 


cys 3S) 


— A r> Uh 4 ) {■ • - 

O 


b 

<>os JO 4- 


(3-4) 


We would like to add that on the top side .of the interval, we have 
& > 0, and on the bottom side we have £ < 0. The skeleton of the pro- 
file with a finite thickness becomes the carrier of the singularities 
(sources and sinks). As a further analysis shows, this is not a pure 
geometric characteristic of the profile (in contrast to an infinitesi- 
mally thin profile, where the skeleton and the profile coincide), but 
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instead it depends on the primary Jf low field. Instead of equation (2.24) 
for a profile with a finite thickness, we will. consider the expression 



where V^ + ' and V^_ are the values of the velocity on the upper and /628 
lower side of the segment. -1 <_ X <_ 1 of the. real axis. Just like equa- 
tion (2.24), the equation (3.5)'isi a requirement that the skeleton lies 
along the stream line. This requirement is only correct for first-order 
terms. In the case,;of the second-order ' terms , deviations have the ef- 
fect that (3-5) must be replaced by the form '('4,25) which is correct 
for second-order terms. This will be discussed in the following chap- 
ter. If we substitute the expansions (3.17); (3-23), and (3.4) in 

♦ 

equation (3.5') s and if we compare coefficients for the same cosine terms 
of the angle £ , we obtain a system of equations just like in .(2.25). 

This system will be discussed in detail in sections 4 and 7, equations 
(4.25) and (7.18). The corrections to the second-order terms are also 
considered there. 


4. Basic Relationships -for Calculating the Plow around Profiles with 
A Finite Thickness. * 

Let us consider the profile with a finite thickness. First, we will 
investigate a profile shape in the image plane Z. Only later on will 
we transfer to plane z, using equation (2.12). I 


From expressions (2.23) and (3.4), we can see that the velocities 
induced by the vortices and the sources are composed of first and second 
order terms. The first order terms can be looked upon as velocities 
which are induced in the plane Z by vortices and sources. They are 
distributed along the real axis along the segment' -1 £ X- £ 1, with 
the following densities 


t 


f ,Y = 2 ( g 0 co'tg i i R- siiiwAj > * 

« ii-i. ^ / 


(4.1) 
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3 —.2 |so cofcg + ff 0 tgi 8- -|- f; & sin 


(4.2) 


In the following 5 the second-order terms are added to the primary 
velocity field (2.17) in order to simplify the notation. We there- 
fore have the task of estbalishing the relationships between the pro- 
file shape and the circulation and 'source distribution along the pro- 
file skeleton. The profile skeleton in the plane Z is identical with 
the segment -1 £ X < 1, Y = 0. The velocity of the inhomogeneous flow 
field into which the profile is placed is given by the expression 


V ~ IV I 4 ’ f) (i^ — oj'cos? nfr! 

H-i) . t ' j 


. 4 - ' 

+ — JS 4 +. 1 . 
5 


Vg v 2 — B z 


~ is H 4 i (kfc A-i s 4 8j| + (I h, -| 

k [>“0 A'^0 V 


4 ■+ * • 4 So 4 “ffij + ^ 7 ii (So - Si,) 4. . • • f ' 


(4.3) 


(4.4) 


\ Jh '■ '••')("• 1 i*). 1 

i ' + i(^‘ + T B ‘ + " ; )( fc + •“.< * 

: K- ft + i [J; Mw +.-v W -8.j] + j (i>( + ; (4 ■ ^ > 

' ’ i > . < 

; +yiirt--'.-].[so 4 -|4 ’Si 4 • • j 

! A 2 “ H-2 + — [ SS' V fc-$2-e 4- SS- ( V fc ,-$*+2 4 v S+2 A)1 4 • 

' J Ll' rt O fc“0 ' 'J I 

l , * ■ . »£'*'..* * i 

*^lrf t*) + f J 

■* . i . . » . i 5 . . n A' , . » 1 j “ 

\\ “ v o “.^0 - (ffb B 0 + 4 + *<+ • • ■) (i + 'j 

j •+%' 4 - [-8* 4-~A 4 : v )(ffo -So/+,^S 2 j.-y^Si-+ • 

: i [J jWb^+ A + x4fe+i fe]+ f ( b 2 4 iii+. ... j> J 

* ^So 4 So 4 “?ij 4 1 8* 4 ii 6 + . s j j|so — So4 “■ "g" Bi'ffi 4 • •• > i 
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In the following, it is advantageous to introduce the complex plane 
£ and the complex function 



(4,6) 

/630 


This function maps the unit circle, in the plane £ into the interval 
-1 < X < 1, I = 0 in the plane Z, 'which is passed through twice. 


In the plane £ , we will introduce polar coordinates p and & (Fig- 
ure 5) : 




Figure 5 


if we also set 

where e is the distance of the point £ from the unit circle j7 Th en from 
equation (4.6) we. obtain 


t £ 


pe* 






»- .. ... U 


/ h O) 


or 


L '•! * 

'• ' ,. i V V 

, ; i+rJ' CM8 \ -:j 

’ ' — i-bin8. ' 

> ’ , l 'S' 1 1 ' ■ , J 


(4.9) 

(4.10) 

(4.11) 


The velocity field in the plane Z, which is composed of the inhomo- 
geneous primary field (4.3) and the velocity field of the vortices 
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and sources, can be written in a relatively simple form, if we map into 
the complex -plane ? using the function (4.6): 






; 4Wo + '2d 57 — 3 - S (2- 4 ~w. 4,4. 

1. \ h t 1- ■ . -5 1 kti* : \ S’? v 

i L i_«. !_ ■- . • ; 




(4.12) 


The velocity is decomposed into a radial part and an azimuth /6 31 

component : 


(4.13) 



On the unit circle (that is, for p ='l), the velocity components 
are given by the following relationships: 


f Oj' ( '* 1 / y ’ ; | 

yj \ — J/rf + ■ V g„ cds-nO -f -- 5) ,4- ; 

,tr 0 ■••,!••• ■ , 3 ■ , 

* Y 1 * * 


' < 1 


( j i f 1 " * * 4 


n=-l 


* 'W* 


; - , +,&.totg la- +. 5i/ n sni )ii> +• 

; V ' * ,«cr0 * 2 1 

' ■* • * ' ; 1 

1 ' ■ ' - 1 + 4 - & - 'S 2:l 'oos 

)i=i J . 


t..; 

> . < 


(4.14) 


(4.15) 


Using this notation, equation (3.5) ’is equivalent to the system of 
equations 


-n- — 1,2; t . i ■ 


(4.16) 


so that the radial velocity component on the unit circle is equal to 
the expression 


! r ’p - slil tf f 2 0 cotg- ~ i) 4- §„ tg ~ $. + £ {/ u Bin jM)] . ! (4.17) 

| t . V ' . 2, _i ^ .«»u • , . / j 11 

which corresponds to sources in the £ plane, which are distributed 
along the unit circle with the density 


ir 

,fj y-gsm# 


(4.18) 
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The profile in the plane £ is given by the closed stream line 
ABMNA (Figure 6). Its shape is derived from the condition that the 
total flow passing through the unit j circular: arc A^^ and the amount of 
incoming flow through the section A^A flows through the section BB^ 
because the arc AB lies along the stream line and the normal velocity 
component is 0 on it. Therefore, we have: 



(4.19) 


i 

i 

i 

i 


V 

■U 

Figure 6 


The value of the first integral according to (4.17) is given by: 




1 


. . jin] the calculation of the other integrals, in equation (-4.19), 
we will use the expression (4.13). Accordingly, the values of these /633 

i 

integrals are given by the imaginary part of the following expression: 

! 

j (4.21) 

(continued) 
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] 1. ,(t?p44Wj d'f>= - j(i ~ S’(m— 





(continued from previous page) 


- ss.iife Hr r' W r- « ir> -W- S& + *H.\ 1 

• ’ ■ ’<“1 t ' . „ ! '1 

• 1 - • i - . i - r. . : 

. !,p-^+ 2 .i l dp, — ' *• ' 1 ' 

= ~-j(l 4 W 4 %) £*•©-» + - 1) 4 ~ fa - fa>|( s * - ' 

: . 2 I < • • ' if < , ^ . Ll ‘ „ '! 

;. + T *0 ^ + f 4 .+ ‘H\ £ ( F« -«■•>[, TTI ( ' (i f-1 

' , ■ " 1 7 ' . * j.* ; ' \ .... ‘I 

-p S )«+*~-i} ^ (( i -I- 4'T 1 1 “ 

i_ ~ ((,i -j- — 2,(^ 0 , 4- ijr/o) [In if. 4 s ) ~ 

, n —1 , ■ J ‘ ' ' ' 

U ((i + s)- 1 - 1) e> & ] - 2& tin (3 + z) -\- ((1 + - 4 <*'4 - 

[*— (i?t+.iffi)fin(i 4 4 4 ~f|i + z )^ •' 

• 1. ' ..4 * # T? ; : ; 4 ; .-i 

!it~ —hi + z)^'~ ij e «»-^4 zK K.;n/ * 


•p_J^hi4 4- 

L «-i . . 


(4.21) 


Since e is a small quantity of first order, if we ignore higher than 
second-order terms, and if we consider equations (4.16), we find: 


) t <> »* ^ ) 4 * f OO 

{ [’ «s4 = s 1 -i —i— si ill) -{r s siii-W £ fl n cos’ 1(9- + 

1 Ji • . A , 2_ 1 4- s i ' , *'” n 

1 ’ * 1 ' »” i 4- C'OS 1) 

j 4 *?»f — S>4 2 (? 0 4'ffoHosH ~ eon*® + 


f f } ! \» t. 


For <?= 0, vre have: 


and for £= ir we have: 


rvfc* 77 . •’ 4 

| 5 ;V^ 

| j 1 < > f ' , 1 « ♦ n # » ' J » r|(J 

^T7 "’/’ Tv 

' C1+«S .t - '4 

},\, ’ ^'dp — -Oi ; 

ijl ■■ -■ - ' - 1 


(4.22) 


(4.23) 
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(4.24) 


For j?=* 7r , the left side of equation (4.19) is equal to 2g o e Q as far 
as (4.20) and (4.13) are concerned,- because according to (4.24) the 
right side is zero. Since the term 2g e is of second order, in our , 
theory we must consider it. Therefore, the equations (4.l6) have beeri- 
given corrections of the following 'type : 
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-* v p ” k Oj : , 

t *“ . F f , * , * ’ * 

* ‘ k t / \ • 

\ 4 (h + 2o<> s o 4 

I' . 1 

■' ’’ >> * 

■ 4 k — 2;- 3,- . . 


(4.25) 


To the integral (4.20), we will add the correction- 

{ A«p <1$ — i/o s„ • ^ {]/ + c'os H silj !)(hi = — g (i e,j (1 — Cos 4) 

.0 ; ■ e ' . . • , / 

■ ' / ■ — jl/d £ o (1 ~ cris^'j-j \ .. ... / 


(4.26) 


so that on the left side of equation (.4.19) we obtain the following 
expression for the total current through the arc A^B-^ and through the 
section AA.^: 


\ = [> ■- 3 . 4 sin 5 > 4 (? 0 02„4 - 2 a j siH ~ 

j -j . ‘ 1 • f ‘ ' /. (4.27) 

■ — ~r 1) 7 s ^"-sin#9 -j- g 0 £ ft (1-fcos^) — — tj n Sp (1— fcos.29-\* 

1 n-B » . n 4 • l • 

By comparing (4.22) and (4.27), we obtain a relationship from which 

we can determine the profile shape; that is, e as a function of the 

angle & . Since we are interested in the' profile shape in the Z plane, 

this relationship is converted from the plane £ to the plane Z using 

equations (4.10) and (4.11). Let us call the coordiina Fes] of the points 

on the profile (X , Y ). Then according to (4.10) and (4.11), we have: 

1 ? ]? 




(4.28) 


After substituting in (4.22) and (4.27), and if we ignore higher /635 


order terms above second order, we find 


iy- j i 4 % i*tf COS ilk + k - go + 2 fa, *4 2d) Col & — fj.k cos 4 

,‘l o , » »« 2 - ' * m ' ,J * 

„ . . ' J 

i+ k ~ + S. 2n sin-^J = - 2d 4 i &) siii a- 4 AT& 4 % + 1 

^ l , * ' , ' S 

4 - g B ] sin' 3» - i S — --4sih«t4.%s (j (i4cos^j~ - ’/ 

^ >7 . ■ . 2- ' n, , • . . / 


"• — .~r k '4(1' — 40 S 29;), . 

‘1 . 

’ | 

X* == — 1 4 —-(Yjj/smii) 2 cos 4. . 


(4.29) 


(4.30) 
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In order to solve equation (4.29), we will write the function Y in 

p 


the form of a trigonometric series: 

Xp = 2i A n sin«a + £ 1 d 2n ; (1 • -&• cos 2«i 4) 4 G^i [cos & 

1 >1-1 di-i l 

’ — cos (2 m 4 1)4] - 


Prom this equation, we obtain e Q according to (4.11): 


, 0 fi-i-lj2M=ii 
l 2 14 s 0 J o 


liui - £ »A a 
0/ siii^ „ti 


and therefore. 


i 

00 -1 

00 

=0 — I>4 tl 4 4 

Ji *-l 


(4.31) 


(4.32) 


(4.33) 


After substituting expressions (4.31) and (4.33) in (4.29), and 
after comparison of coefficients of the same sine and cosine terms of 
the angle ^ , we obtain two systems of equations : 




V-'. _ : ’ 

- do do + -J di 


1 ^ \ 4- A 'do ~ % + H + i-Ji - 4 fra] 4,4 

\ . ‘35 J 1 


. .A 

. '• 4--” 4. 2 (2(^ 4, 2q n 4 !h . 4, tf 3 ^ 4 :j)'.4 ' ' v ; 1 '.; 

, , • ‘2 • “ ■ • • -V ' >•' > - .-i 


1 ! 


4 




9 2 A„,(~ ffrf-1 + 4.-1 + dn+1 4i+tj 4- 4uf ^ £ &a»i4 £. C^m+l] +, 

■“ «-* 8 L #<««l ' * ni»»i £ . J’ ■ 

2. Oii + “<?.]+ 7^4 2, t? 8m+ i+ o 3 ]- v £ ^'[0,,-i - buijj' 

m-1 ^ J L.w.,1 'J ^ -,ti **3 

* * r 1 ’ , • * 


, ' 4 do 


i 1 ~ . . 

+ ~ £ (“2«--^4;4rf--2'-“r %4z — 4Hi4+ , ‘ , < 

« ««J " ' ' > 

<x> “1 ' 1 i * , • t 

^2 + 2 2 °4 +ifc- 2 /^ 2 rti + (*7x H^) £ ' C4<+i 1 

jtrt. — T 2 1 


^ ]|«l) ’ * * * i * 

■ i ' . 1 •; , a , . 

^2 di) — A n 4 -^Ai ( —q 2 -|- q A 4 — n 4 ) -4-“, 

.+ ~^A 2 [ 2q 0 -j-22ri + db + ~4o 4 •~dn 4 4 — ?4 -** ,>.• , 
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* 1 i 1 j| ^ ^ 

. “'2li + -±- do + — da — A 2 4 — A\ (2'g,; -|-'2i/ 0 4 4j !. "4- $ 3 ' ~ di) 4 1 

2 2. 4- 2 7“ - 

v * 2 1 v* ‘ i *» 1 ' * ,' p k 

+ ^aUo— §a + £o 4 -j *?4 774 M 4 ^“^^0 +?r 4 <H — ?s^ 4 - \ (4.34) 


l( continued) 
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(continued from previous page) 


v-l,'’ ' ■» , ‘ 1 , , • ?t< ' '* : < .. .1 • ‘ v Ai< j (5 

>\ ' ' • • "4 “-^4 ( 2 2<j + 2 ?o + 2 t + &. — -K)' “i" ( 

* £ * v r * J | 


’ ct ^ J ^ tl ^ 2n~3 ,t^n~3 ^tt+3) *f _ ffo 

^ tt~G . * . 


. n -i 

+ \ (On. + (fi) S 0 2m+1 - i (r/j -f g s ) C. A - 

2 'm-= J. 2 

11“ 

~tOg t-s 4- — y (f/ 11-3 ■ 2» i- 3 ) 


«?a+ V «S °* 


, t 

’ 1 


2 •' 


o 

u 


X 


X 




i(i 4- 2u~ ffu+;W);§^ 2m + (:&+ & + v4 S ^1+^- £ ({?» ~ ft,) & 4 j 

!■ t ' ■ “> • .rt*i • \ -2 ; w „i , r , (( -S-rtTB' 1 ‘ 

V.X '4.2 0 y.Xh f4?4" .ti, On — •rpffb 'fj, * 4* f ■ r ‘ v • '/‘•if-’.l 

(1 -yq 0 - q 0 -f uj} y C zm+l + .(2q 0 + 2 % + fe)f y i % n - ~ (J 2 ) ^ " } 

.. Irt»l . ‘ W-i 1 i . . 


•■jI 4 
’ . 2 


(Sa X ^ 53 1 ^sw+i ~ d ;i j 4- 




•v +t£ ~ Ki) ^n+i) + ffo [^L a -Jr s 5 

2 n-3 V ?»-a J : 


"t y {-^n 2n4-I "t ~,2<) y 1!-i n } r 


Vv 


0 


tt-1 


•- U + 2o - &• +.%) b t + (&+ 2d + **j ( | Gy + t - <4j^ * | 

-ia»-Ps).f s 5 bwHa - ^t'. 4 , ■ ■ -I 

- \»-i / 

* * , . ► I 

~rS (2« hj + d\ t+ 2) ,-r i/ojfx ^ 4- 2 y ^ ii/i'4- j 

2 ».~4 • V i - «=.8 / . 2 ; 


t" ' y (^n ffj »+2 "f" r“ Oo y }i, .-4nV. 

2 tl „i , 4-' „_i ' ‘ 

■ ■ ' ;* >• i t' • ' * • ! 

r — (4- 4-26 — 2li + ^0) ^3 “ 4- 2b -1- -j^)x' 2 -H ^4) “ ■ ’ j 

i ** p * • k bo ' 1 * \ , 1 f j, , \ ^ J ^ 

' 4 .'■ ’ o’ .,{* 2 ) (y ^s] ,( 2 s " J ",i A d) I yX 2 ' , *~~‘ tp ^ 0 ) ■ .4 

• _ "2 • - ; •• • '\£3i ■ ■ • ; .%• .)■•’, *.4 , ,h 

*’ * I ' * ' , ^ * 1 ‘ ' t I I 

: -v(r/4-^) f y c/ a#l+1 + (2„- &,);<<?, .-*4 C?; +a J + ' ; 

•' 'i 2 /■ \.>»-»i ’ J, 2 ji-s ' • • i j 

, J - t ' j 

’+ 2b ( 4 . 3 + 2 y 4„| y J, a 5f 3 _„ + y tA,'ffrt-,-3 4 - On) — 0 , < 

V h-4 / 2 «-.{ £ u„, j 


( 4 . 34 ) 
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( 4 . 35 ) 
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It is often advantageous to replace the first equations in the sys- 
tem (4.3*0 hy the following equations: 

‘ ! — £' ( 1 + 3? 0 + — Y> + Yi t A « } * 

2g 0 - £ 1 - ([„ ~ 3§t 0 “ | {— i) n 2« f £ (-IT ?«.]• ' 


These equations are obtained by multiplying each equation of the 
system (4.34) with its own order number. Then all equations are added^g^g 
once with the same sign and the second time they are added with the 
alternating symbol. 

The first two equations of system (4.35) are linear combinations of 
remainder equations and therefore will not be considered in the fol- 
lowing. It also follows from (4.35) that the coefficients C n are second- 
order , so that in systems (4.34) and (4.35), all terms in which coef- 
ficients C n are multiplied by any first-order term can be considered 
of third order, and can be disregarded. This means that the C terms 
drop out of the system (4.34), and the system (4.35) has a purely dia- 
gonal character. 


(4.36) 

(4.37) 


5. Profile Shape, in the Vicinity of the Leading Edge and Trailing 
Edge 


As already mentioned in the previous chapter, the Y coordinates of 
the profile are given by the following expression: 


3T„ — £ A n sin-»9 + . £ C 2nt [i — cos 2 •>«$] + 

- ’ * i i 

+ £ '‘C f 4 m+ x [cos& — cos (2m + — tcsS&jS; 7t . ! 


(5.1) 




According to (4.10), we have the following for the X-coordinates 
of the profile 

i'j. -fl. + -i{Y J1 /s!n^] 1 co8d, - (5.2) 

. . ‘ * L J . 

or after substituting equation (5.1) and ignoring terms higher than 
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(5.3) 


second order: 


, j,i = — CtiS ^ I 1 -j — -p- JLj -f - dj “H — r d s “I" dj d 3 ~|~ iL^dg -j- 2d,; d,| -{7 ■ • tI” - • 

+ ^Uid 2 4" -dj d,j: -j-*Sd 2 d 3 • • -3 I‘o§.$ -t - [do-f-^dj -jf 2dj^ dg^-ti^-g d 4 -j- 

4-2djd- 3 -|- . • .] eos2& -j- 2 [dj d 4 -j-d a d 3 + : . .]'cos3 5> 4- [d 3 + 2 d a dj 4- , 


4 : 2d , i.d 5 '+ .'..'jfcoiU'O 


s 


s 


Prom expressions (5.-1) and (5.2), we then obtain some important in' 
formation regarding the profile shape in the vicinity of the leading 
edge and the trailing edge. 


For the leading edge, that is, for # = 0, we obtain: 



or according to (4.36), 


(5.4) 
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% = 0>. ±1 ~ '2qih . ! . (5.5) 


.Without any difficulty, we also obtain’ the following: 



- Bz _ 


and the curvature radius of the leading edge 



Expression (5.7) can be written in the following way: 


•%== fff- i’4- H + i - s & + ; 

\ “ K * . X-2 , tt-o * 1 


(5.6) 


(5.7) 


(5.8) 


Finally, we can prove that for 

X p = -1, we. have Y v == ± ± itj. . (5-9) 

From equations (5-5) a (5-7), and (5-9), we find that the profile 
shape can be approximated by a parabola in the vicinity of the leading 
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edge: 


‘ jjT=’± + i 4- 2ri) 1/2 = ±,^(24,+ 2' +.R i )^:i (5.10) 


Similar conditions also hold for the rounded-off trailing edge, that 

' % * 0: 

' JT=o; i,-t Hp; ( 5 . in 

, \ ^ L #«1 J . , J 


is, for #= tt , where q Q 0: 


or, according -to (4.37), 


Also, we have 


!;-<>; ± v - 1. + 225 • ! 


aiy . 

f - OO . . 


iix. 


The radius of curvature of the trailing edge is 




or 


! 


k co 

IS 

Tt* •• 

( —I)" n dj 


f Jf, 2 — - 

. 1 + — 

■ 2 

f r (-l-)-*Xf 


* , i 

i— <f.i - ( -!?«*.+•&( -it FifJ 

. j • ' »-s- - ' • tfMo ■ , JM 


Also we have: 


-i, n = ± -*55 = ,± K , } 


(5-12) 


(5-13) 


(5.14) 


(5.15) 


(5.16) 


so that the profile shape in the vicinity of the trailing edge is ap- 
proximated by a parabola: 

j i 9 - ± ^ - ■*;)“ -j ( 5 - 17 ) 


Figure 7 
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In the case of a sharp trailing edge, that is, for # = it and q Q = 0, 


we have: 


and 


X„ = i,'r r — 0 

d f 


dX„ 


=0. 


( 5 . 18 ) 

(5.19) 


Finally, we will make the transition from the image plane Z to the 
physical plane z (Figure 7): 

- *■ i T 

z v “ + * 17 (^p) — / (2V)]. . (5.20) 

The square bracket is replaced by a Taylor series: 


r }(x p + iY v ) -J(X V ) + ... ^ 4 


cli»’ 


(5.21) 


and since Y p as well as dy g /dx are small, of first order, we can re- 
strict ourselves to the first term of expression (5.21). Equation 
(5-20) therefore has the form: 


<V— «» 4 i (I/p ™ Z/») «*■ < "V. (i + < j' 


In addition to (5.22), according to (2.12), we have: 

* =* XV i 4 •— (Xy/sih '•&)*'] cos ft. , 

• L - * - J - 


(5.22) 


(5.23) 


If we decompose (5.22) into a real and an imaginary part, we find: 

Xi, - : i , - X„, ■ (5.24) 

els? 

which means that the thickness t in the. plane z is plotted along the 
normal to the skeleton, where 


t = r„ 


1 4 


r dj/« V 

l da; . . 


(5.25) 


Since the root value in equation (5.25) differs from 1 only by a 
second-order quantity, we can write the following for this theory 

( 5 . 26 ) 

Therefore, it follows that for the profile shape in the vicinity of 
the leading edge and the trailing edge is the same in the physical plane 
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z as in the image plane Z. 


6. Velocity Distribution along the Profile Contour 

When calculating the contour velocity, we will start with the velo- 
city field in the plane £: this velocity field is given by (4.12): 

“s - \ (i - £) J 1 + £ < ■ - !>• (P. - ”,) | «* +E-) + 

. + 2( ?fl + <jj„) IX + I)-' + -lt h (?-!)■' - £ (j, -I- )\< f,) (-?)-” 

Prom (6.1), it follows that the contour velocity in the plane Z is 
given by the following if we ’dividei this expression by the derivative 
of the mapping function (4.6): 


tl2 



( 6 . 2 ) 


and if we substitute the following for the variable 

i; £ )h- ; » ’ ‘ (6.3) 

where e is given by relationships (4.19), (4.28), and (4.29) /642 


JJ K i Y k _ ] + £ (|x n — il n ) ~ [(1 -f e) n e in5 + 

2 




+ (1 + e )-»e<"»]-( 2o -f^)* 

2 . * .3- 

ffo 


(6.4) 


(1 + E ) ©-» - 1 0 (3 + s) o-« + 1 - 

~ S (ft. + ^Tfc) (1 + «T* cM 


1 


If, in equation (6.4), we substitute v n according to (4.25) and if 
we ignore terms higher than second order, we find: 


oc 

U K — iV K = 1 + 5j ?» ( cos ^ «$) - i g 0 (1 — sQ -f 




+ ig 0 s o cos & — (go 4- ig 0 ) 

(3 4- s) (cos il -j- i sin &■) 1 (1 -j- s) (cog & j <uh 4) -f- 1 


j 

(1 + e)(l + COSS) + 4 s 2 
2 


2o* 


+■ 


(1 + e) (J 4- cos£) + — s e 

2 


(6.5) 


(continued) 
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|4 2 {#„ -j~ g„) (1 — s)' (c6s 9 4 i-siii kj — - b e), (dos. « 9 -M alii 9ifj’4 j (6.5) 

4 • ' a-s . , . 'x 

. ‘ g„ (sih49'4- iii-e 'cbg ! n9)v- t > \ •>. . 

. . ■ Si~i , i . i ' : , . * <-4 ./ 


By decomposing this expression Into a real and 'imaginary part, we 
obtain the two components of the contour velocity: 


• ‘ • ttk — 1 + 53 ?„ COS. 9 9 


Ho 


‘.(1 4 s) COS, 9 — 1 


Ho 


(1 4- 5 ) cos -9- 4 - 1 


■ ;/ 

“hr* I 


[■'" (1 4» (1 - cop 4 4 - — s 3, (i 4- e); s {‘l -!- eds/9): 4 4 .s 2 ' ; 


r .. « r , 

f ’ •- -k -i- tfb)..(a - e), cbs 9 £ 5„’(i — iuj fats ?i9 4- . 

y; ‘ ’ n -i - . ■■ ; . 


k" • 4 Ho t 


si h 4 


1 I , / 

7 


1. 

k’< v ' 


i pi-.* .r', vf'"4i 


( 6 . 6 ) 
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^4— flbr 


.e 4 -T s 2 . 
2 


' it !/0: j 7 ( lo r 

V , ; . ;'(i 4 (1 -r-iios'4) 4; 4 * A • . i 

■ 2 t • 


siii ■& 


j 4 » * , * K 


cOg. 9,4’ 4- ; e 2 
v ■ 


.4 


; 4-3b — 


siii' -9 


t ' 
\ 


i -tt ‘ cbs 9* -f- — 


2 4- $j) siii 9 4 £ 3b sin «9 




— 3o s o (i 4- cos 9)' 4“.'s (41 + —■ 'So 4* 2 § 0 ) sin 9 4 


• *, 


■ '4- e % M?,, — i i) jsin 7ik E £ 4; g H cos *4. . 

^ * * * ’H«-i « * 




The contour velocity is then 


(6.7) 


( 6 - 8 ) 


Beyond the surroundings of the leading edge (and for q Q 7 0 also 
beyond the rounded-off trailing edge) the' contour velocity differs only 
slightly from 1. Therefore, -in this region we can derive a simple 
expansion for the contour velocity.;- Prom expressions (6.6) and (6.7), 
we obtain the following by ignoring terms higher than second-order: 


<x> ' ' ’ J 


47* =14 It 4 f k ~ Uo + 3 (3b + ( lo) cos, 9 >- X'. ^-^cds ,n9 + / 

i , * * ff'-U . > “ ' \ t c* '• 


<v> * . 




*4 f fffe (Jdtg‘4 ^'4* ti 3b' siii 494-4 


IH«1 


, .1 — cos -9’ • •• 


(6.9) 


“1 T 


:3ik 


1 4 cos O' 


• 2 (ff„ 4 3n) cos 94 Y, '•Hu cos 9:9 


( continued) 
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(continued from previous page) 


V K = q u c ofcg ~ a + <Zi, tg ~ - 2 (4o + !ZJ si” + 


(6.9) 


4- S </„ sin n a- - f/„ e u (1 4- <•<« ») + 


11-2 


4- E 


J » (P« ~ (In) sin u tf + . - £ i)jf n cos n ft 

«:»! c J — C*OS tr « - 1 


( 6 . 10 ) 


In the following the letter <p is the angle between the tangent to 
the profile (Figure 8) and the X-axis; and then we have 

Wji - tr k cos 0 4- V K sftt *. j (6.11) 


Here, tg cp is the directional factor of the profile tangent which /644 
can be written as follows to the first approximation according to 
(4.31) and (4.34): 


ri y 4 1 * v .-i 

% ? = -r~f = ~ Y> nA « cotg^-i) + rt-S (-l) h 7U„tg— D - 

> Q-JL V & u**i Jj a i 2 , 


£ 2 m A> inl sin 8 — 2 

m«*l X -* 


5] (2 ??t + 1c 3- 1) 4„ )n+/ . hl 

m- 0 


. siu 7ri), 


1 V JL - i 

tg V « ® y = 2o cotg — a- 4- & tg - a - 2 ( ffo + g 0 ) sin {> 4- 

co 

4- £ q„ sin » a. 


(6.12) 


(6.13) 


Figure 8 



dx. 


Therefore, can be considered as a first-order variable, and 

. dX„ 

the contour velocity can then be written in the form 


tl 7 . 


U, 


T—i 

i 

f dX„ \ a 

C4 

1 

Idx ; ) 


. dX p 
K rlX. 


(6.14) 


After substitution and if we only consider second-order terms, we 
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find: 


]V K I 4- £ \J.„ rox i>8 + g„ t - r/ 0 4-2 (</„ g 0 ) cos 8 - f] q n cos 7)8 -|- 


1 


4- ih -^- 8 4- £ 9« siu '»8 + £ 

_ 00 * 1 td Y 

~ ^ (5fn 4“ ffo‘) 40S -5)- 4" <’OK 9} ft 4-, ■■ ■ ■ ■ p 

st»a J IdA*, 


£o j So 


1 — Cos 8 14- 80S ft 


So oofcgi-S-'-K (6.15) 

' Ji 


4- fij t.g -i- ft - 2 (s« 4- So) ft- 4- J n-g tl sin nft - — f ---P-] , 


ar. 


where the value of the derivative dZj, is- given by (6.12) or (6.13) /645 

and e is given by: (according to 4 . 3_li.il 


S — ^ -Agm-f-i 4 - ^ 


It .1 


^ -42)11 > 11 -hi 


7)1 *»0 


COS 11 ft, 


(6.16) 


or 


(sc 


4b - q'o 4- — ?2 4- - (Li + . . • -I- So 4 So -I- T flfa + 


3 


15 


44 


J 


•1 


+ 4 ?6 4" • • • 
12 


12 


)“*-(? 

4 


Sa- 1 Sa + .v.| cos 2 ft 


T 4a ~ ~ & 4- • * «’l cos 3 ft ~ Si | 


cos 4 ft 


(6.17) 


( i , 

— i ~ ?5 4" • - • cos 5 ft 4" * * 


• ar. 

For the derivative ^ and for £ we will always have two expressions. 
When solving the even problem, we will use expressions (6.12) and (6.16). 
The contour velocity is then given by a linear expression in terms of 
the coefficients g n and q , which is important for calculating the con- 
tour velocity for different incident flow directions. When solving 
the odd problem, we start with the selected circulation distribution 
(coefficient g n ) and source distribution (coefficient q ) . This selec- 
tion is specified on the contour by requiring a suitable velocity dis- 
tribution. Therefore, there must be a way of evaluating this velocity 
distribution before carrying out the entire calculation. Relationship 
(6.15) is used for this, in which we substitute expressions (6.13) and 

(6.17) for dY /dX and'e. 
v p 

As already mentioned, the expression (6.15) for calculating the con- 
tour velocity cannot be used in the vicinity of the rounded-off trailing 
edge. There, we must use formulas (6.6) and (6.8) for obtaining the 
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contour velocity. Directly at the leading edge, that is, for £ = 0, 
we have the following according to (6.6) and (6.7). 

Uk - 1+ £ K - H (1 - J + + Q« U ~ ~ *ol - 

"o U • / l 2 j (6.18) 

— £ tf»(i ~ >»«■)* ■ 

; * 

• * - .{ 2 \ 00 

* _ He — //ft + 1 ~ 2 toj — % J] ngy, >. < (6.19) 

If in equations (6.18) and (6.19) we substitute according to (4.33) /646 
we find 


and we have 
• Vi r 


If in these expressions, we substitute for ZnA n according to equa- 
tion (4.36), then we have the following: 

u * “ ~ (77 & + j & “ S > 

a- = -i 1 + 32 b + «o — »to + Yi V-h ) —-?b £ *to - 3 2 0 (fo . 

It should be realized that terms higher than second order were not 
considered in expression (4.36) for EnA^, and that EnA n itself is a 
first-order quantity. The ignored higher-order terms have an effect 
on the reciprocal 1/EnA^ already for the first order terms, and there- 
fore expressions ( 6 . 23 ) and (6.24) are only correct up to first order 
terms. In the case of a flow without circulation, both components are 
small of second order in the case of smooth entry (g = 0). For a 
non-smooth entry, U K is small of second order and V K is of order 1. 


(6.23) 

(6.24) 


£ ° “ £ nK , — 


• _ 1 ±,i® 

» O ^ | L }iA ii 


+ 3 + S ?» + 3 «o 4- — S a.. — S n A » — to + 

53 n -° n = 2 "-- 1 \ 2 

3 _ ' 06 

' + -T to “ 13 ’ 

& J 

Vk = to j -i £, »*•) ~ £ £% fc • ‘ 

53 1lA n 2 - 1 ' *-» >=-' 

\n-l , , / . 


30 



In formula (6.12), which can be considered as a formula for the 
contour velocity at the leading edge, it is often advantageous to 
replace EnA n by the radius of curvature of the leading edge, accor- 
ding to (5-7) : 


Wj, = V K ~ 2(j n (lii 1 _ ) - :Rn 5 «</„ . 

1 


(6.25) 


In a similar manner, the following expression gives the velocity 
along the rounded-off trailing edge (6.26). 

, , 1 

£• £ "4* (6 - 26) 

• ' \ ir-d 1 7 


The X-component of the contour velocity on the rounded trailing edge 
is small of an order higher than second order - , and therefore can be 
ignored within our theory. It should only be realized that formula 
(6.26) has a theoretical meaning. This is because in a real flow there 
is a boundary layer separation in the vicinity of the rounded trailing 
edge . 
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The formula Is derived .for the distribution of the contour velocity 
applledin the image plane Z; therefore, we must derive the relationships 
for transforming values from the plane Z to the plane z. For the velo- 
city in the z plane, we have 

l_J “* (6.27) 


w 


where 


dZ 


a z 
dZ 


TF } 


is the absolute magnitude of the derivative of the mapping 


functions (2.10) 


<lz 

dZ 




a/ 

d Z 


( 6 . 28 ) 


On the abscissa -1 <_ X <_ 1, Y = 0, we have, according to equation 
(2.15); 


a/ • 

dZ dx ' 


(6.29) 


and since the profile thickness is small, the derivative df/dZ of the 
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mapping function can be expanded into a. Taylor series on the profile:' 

' % = ?*+ ! (6 * 30) 

: dZ dX y das : 

If we only restrict ourselves to terms up to second order, then we 
can omit the difference X__- 1 and all other terms not shown of higher 


order. The derivative 

dL 
Ax 


p ;2l 

Ua/j 


Oilcan be written as: 


f ~ A 


V <1 £cj) d> 

1(1 coj 


rU) " sin ub ] 

-r— — — V nB a — - — ■ ! 
tlx ■ < sm-ff ; 


(6.31) 


After replacement in equations (6.28) and (6.30), we find: 


:||>i *t’s<M** • + jj * * 

‘ClZ r ' r'.h««d« 1 * ri—i i ' SfiSw * r/‘ 


iiix n& > 


(6.32) 


For an accuracy up to second order, .we therefore substitute the 
following in equation (6.27): 

• 1 m = i -■ | ’4 ~ i Ik +■ . ;) - M +| 


: jm 
h t 


w 

i 


; + XX + • • b - |1 k + XX + | •] -j 

‘ ' J ' ! ; ’ ’ . . ’ ' ’ ' ' • * ' ' < 

-(x H, k X -I- | Xk -T- • •] <*s to - k 4 B a ty 4- i 


+ ~ % X + - » w k - .] B 0 & 5 -s- -J % 4j + 


I ' 




i 

1(6.33) 


2 ' 


i ■+■ — i? 2 i? 3 -b „m.| bbhM — 4- M. s h 2 - + . i .) sii! -& — i i 

; • , 3 ; ‘ } ■; 1 .■ ,Ju . y . ' ' '] 

— (Mi B 2 + J.-8 + &A* k +'2l, X + - - -J'sim 2-a — (3/i x k +■. . / 

+ 2lilh + -ib Hi .4-- * 0 siii 3«. —'(Lii Hi H- 3d 2 H, + 

: ' . . + Ai X + s . . ) silt k + / 


v± . 


The contour velocity is then given by the product of expressions 
(6.8) or (6.15)' and (6.33). 


ms- 
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-SECOND PART 


FORMULAS EOR NUMERICAL CALCULATIONS- AND EXAMPLES 
7 . Formulas for Numerical Calculation 


In the previous chapters', we used trigonometric expansions with an 
infinite number of terms for the variables to include the general case. 
In practice, it is sufficient to use relatively low-order 'trigonometric 
polynomials.- In this section, we will .discuss formulas for practical 
applications, in which the first few coefficients of trigonometric ex- 
pansions are considered in the first order .terms (usually-, it’ is not 
necessary to take more than six terms). In the second-order terms, we 
will only consider those terms where the' index sum is smaller or equal 
to 4 . 


/649 

Geometric Variables 

* 

- The coordinate system is located so that the skeleton chord lies in 
the x-axis, and the origin of the coordinate system is located at the 
chord center. The length of the skeleton chord is called c*. Instead 
of the variable x, we will introduce a trigonometric variable & with 
the following equation on the skeleton chord. 

' '£C" — — 4- COS -ft; jj -<!'{)•< TC 1 (7-1) 

Li, .. LL.,:. .7.. “Jj 


The skeleton shape is given by the trigonometric expansion of its 


derivative 



(7.2) 


*In this chapter, intended for practical applications, all of the quan- 
tities are dimensional. The skeleton chord length Is the fundamental 
length for length measurements, and the velocity U 0 is the fundamental 
quantity for the velocities. In order to simplify the notation, we will 
use the same notation for these dimensioned variables as in the complex 
plane z. Since the theory is not required for performing numerical 
calculations, there is no danger of confusion. 
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where 


t 



From this expansion, we find: 


where 



(7-3) 


(7.4) 


is given by the expansion: ( 

_ 


The thickness distribution is plotted along the skeleton normal, and 

rr: i 

(7.5) • 


■ t = - '{ £ A n Hiii -f- y; ,{c, m ,(i rr- co« smU) + 

-I- 0 2hl - |-r(cof? — COS (2m 


The positive values are measured, .upwards and the negative ones are /650 
measured downwards from, the skeleton. The base point of the normal to 
the skeleton has an x-coordinate given by the relationship 


or 



(7.6) 


(7.7) 


Figure 9 
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The curvature radius of the leading edge is: 


“ T IE vA «) ' . 

^ \n~>l / 

The curvature radius of the trailing edge is: 

n ( 31 \£ 


(7.8) 


(7.9) 


The distance between the profile nose from the end of the skeleton 
onthe leading edge is 


\ — — Jt , . 

■> 


(7.10) 


and on the trailing edge. 


A 2 — — ^2 • 


(7.H) 


Aerodynamic Variables. 

The components of the primary velocity field on the point where the 
skeleton is located are 1 given by the trigonometric polynomials: 

: u = U n f .1 -j- £ cos «•!)], v — ; Z7 0 £ % o.Os wlb | (7-12) 

. * , V < n-'-O / } • > -n~ 0 , 


The circulation distribution and source distribution on the skeleton 
is given by the following trigonometric polynomials: 


- y =» f k cotgi ^ + S to sin 


(7.13) 


? = 2tf 0 


So cotg i -9 + ft, tg i 9 - 2 (q n +■ &) sin '9 4- 


J /-1 


+ £ sin ra9 


(7.14) 
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The relationship between the geometric variables and the aerodynamic 
variables is given by a system of linear inhomogeneous equations be- 
tween the coefficients of the geometric variables A , B n , and C , and 
the coefficients of the aerodynamic variables g , q n , p , and v . In 
order to solve this, it is advantageous to decompose the system into 
three groups of equations and to then solve one group after the other: 
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The First Group: 

s 2o " So 4- — 2 1 “ A + jg# ™>2,j 4- — 2 2 4- ho ~ h»j + ■ 

’ ' s> .. + — A z {2(ln 42g 0 .4 ,hi)r ■ 

,i , _ 2, ” ” ■ ■ ; 

'■ j ' i 

X X 1 f *1 ' i * 

p 7o + — 2o 4 ~ 'Sa — -4 2 4* — (2g 0 - + 2p 0 4- <f» 4-hf ~ &t) 4 .! 

i .. * 

> +42 (2o — 2d 4 t x o) + -4t(2</ n 4 2<Zo 4 hi), 

• 1 , i i ’ . < (7-15) 

4 (® 2 2d) -4, — A\ (g 2 54) 4 •— 4h> 4 4) 4 

i> Z Z i 

• ' ’ + 2i?o 4" hi) 4 -d-a (?o — 2n + hn)» - : 

1 ■ ' j . X 

— TT (Ss- tfc) =4“ 4 4j (' 2 * ~ it,) — ~ 4 (2a — 4)’ 4 

t i* -l . a • , 

. •+ ~ 4- ( 2 2d'-+'2® 4 h,4+ 4 (2o' " 4n + hi), ' *; 

L ■' • P ^ fy+iKp S;- ; '•}’ 

It is better to use the following equations, instead of the first of 
the two equations given above: 


•; 2 2i == S 4 S? rt 4 g„ - £ 44 £\hi 

I . l ,i«2 tl“0. " ,) 


(7.16) 


22a 4 S4-l>“+i4, p 2o -8&- £ X -!)•& + £ ( — i>" 4nl ‘ 
.... **”4" ' - "~ 2 ’ »-o J , 


4}.‘j (T. 


The Second Group: 


,y 0 = __ u o 41 v n 4 y 0 g »l n 4 b 0 , {id - - — 4 . j* y. a 4 , - 1 

’ • , , • ‘ , ri"i • 2 2 ■' , ■ ■ 

r &) ;?.f & & + 2 0 ) 4 ^ a (23 p&jr „ .;1 

; 4 a* '! - /• ' t • ; ? ■* 1 v ! 


{h - 24 --‘v- fa fc S *X* ; + i?i fiv+ I it*] 4 i JiaiV- , 

i v ^ I 2 


5 

^3 (ffo — 2o) — lli (2it 4 2o)> 

• 5' 


i (7.18) 


« * ’X' " f c j 

2s = - V 2 4 4o h 2 + 4 4a) 4 7+ h„. — -J {{ft — 2a), , 

~ 4 • o 


...1 

' * .. i 
■ • \ 


< v 3 + Dd +: 4 “-jife 4 — '4a Uf4 i? :i h 0 ^ ■“ ) » 

r ’ •' h , . i 

j , , . 4 . ; * , K 4 ' . , * 

i 2d = 4.i — v 4 .4 £j h< 4 — +1 ha 4 4- R, h. 4 *- 1I3 hi 4 ^4 ho j ; 

! ■ *< ^5 >2' 

1 2n ~ v ii , ‘f! 5^ 5. 
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The Third Group: 


a n -A n 4" ffo f-^2 4 - 2 — -j- {fa QL — '4-a) 4* 

* »•«! \ n«»3 ( / 3 

1 pi = l/() (i-3 + Si 4 ). T- ) i fli.jlj — ~ (J 2 A\ y_ ' • ■ * . ‘ 

r ■ j 3' 3 ' ' . ‘ N 

3 . ’ 1 J ' , . * ' 1 ’ 

h ft “ "£■ ( Ji ^2 ~ ~f v ... . ,'( 


(7.19) 


All three groups have the largest terms in the main diagonal, and /6 53 
therefore can be solved using a method of stepwise approximation very 
easily. In most cases, the first two approximations ' will be suffici- 
ent . 


THE EVEN PROBLEM (THE SECOND MAIN PROBLEM)' 0F PROFILE THEORY 

When solving the even problem one determines the circulation dis- 
tribution and source distribution (that -is, the coefficients g n and q ) 
for a profile immersed into a given primary field. Prom the calculated 
circulation distribution, one determines the total circulation, the 
lift, and the velocity distribution along the profile contour. 


Since the profile skeleton does not have a purely geometric profile 
characteristic (it also depends on the primary field), one starts at 
the profile center line when solving the even problem. This is because 
most profiles are produced by drawing a suitable symmetric profile along 
a corresponding center line. 


The beginning of the center line is at a distanceAj^ == 3'ig,from the 

2 

leading edge (Figure 10), and its end in the trailing edge, or at a 
i 

distance A a = — H 2 for a rounded trailing edge. We will determine the 
. 2 - 

coefficients for the trigonometric expansion of the slope of the 
middle line : 


■da 1 


£ cos iHh, ] 

ii — O * ‘ 


(7.20) 
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Figure 10 


Also, we determined the coefficients A n of the symmetric thickness 
distribution : 

■t — ■> — y a> — — 1 jeosJh r w 

; -3 bUi • (7-21) 

Considering that the nonsymmetric correction of the thickness dis- /65^ 
tribution (the coefficients C n in (7. 5). are of second order) we can 
take the following curve when calculating the skeleton: 

,M»~ ?/m{ [$ 2(1 cos 2& ) -j- G z (co&8 — cox 3-8.) -fa 0$ (1 — cos 48 , I (7*22) 

The coefficient values which are substituted in system (7.18) are 

then ...... 

B'o — B 0 — 2C$ , 

, 4 - 'b;- 4(^ + 2^), 

f Bj = B a T-0Caj- - ( 7 . 23 ) 

. $ - ; ic 4 , 

where the coefficients C are given by expressions (7.19). By solving 
(7.15) and ( 7 . 18 ), we obtain the values of the coefficients g n and q . 

The total circulation is: 

F » Ttcfro • C7.24) 


The velocity distribution on the profile contour is given by the 
following relationship in the region outside of the rounded trailing 
edge (and also outside of this trailing edge, when the trailing edge 
is rounded off) : 
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& 


T^T ; {‘'+Xft' 


5k_- 

id 

l^i 

i •' 

* +ffd cotg — 4> + £ 0 „ si fiH-i) - 1 - £ 

‘ n l 

at 


t. 

cos nil + <Zo— v/V-l- 2 (g d 4"2d)c<i^>— Y, ‘i, cosii# 4-' ’’ 

So I 


iU 


\ 

So' 


1 — COR ft 1 -f COR !) 


2 (?,( +So)^'S^t: i 


+ 5j Jiff,, COR 1(9 
! r n i 


* dx 


So < Otg -^-4> + So tgx* 2 (So' ■+ So) Kill 9 


( 7 . 25 ) 


’ 4 ' ' , j /ly Ji 

4- V a sill ti$ ' i’i 

i ' ‘ 

• ■ ' . -A ” , 2 (IX. 

j . - 


fi T #■ +V\+.|-‘ V | Vi 7J 8 + S« |i>4 4- i Jh'+jJU + + .' i . ] 

ply. / ; " ‘ ’ f ' : 


’ ■• • '■ ‘ . • J*;,:® 
V i\V -*’1 

SXJl 1 


7 .(If = [li + v j7?i 4 Vj |i? 0 + ~ B t v,,/i, -s- So ( Y «* + .£,) + , 

H ■ ■ * * » , * ( * ' 

.. ' i. l ' r» * * 

r /]Ia— ii*i+Vo^2+ ~ V 1 J*i + H & 0 + ~ V 3®1 - 

h r - ,_ / ‘ , ( *‘ 1 

* > * 4 ' - 

* 'w«t»8+ V.A+-J- ^1*,+“^+ ^K+r^ • • ’ • 

2 2 5 ‘ i 

;']t* < -W.+ v tf ^+|-v l JJ s + iv 8 B t '+i^i+v^ if . . ■ 

* »^ 5 ; • -. .' 

| s = A{ -}- X ;j -)- 2 (yl 2 +X) b6s tf’X ^ X a 'co ! s -2& X 2 ^ 4,008 39--; 

; it* = ( coi g — 9- IgA £-2 (2X 2 + Ji 4 ) hin 9 - Gvi „ sin 2,9 - /■ 

> flX„ \2c/ 2 ' 2c/ 2 ' 

■•. ‘i ‘ , ’ 

( ,•■ : — <$X, sin 39 ' . •. - , 

— I = i + i Bl + 4- {Jif +BD if B n Jii + — Bi J? 2 1 (JOS y -I- 4- i 
' 2 4 - \ 2 ) • V 4 

‘ . i4 « ; / . v •;* .* . t , , * 

+ i. i?i^ 3 j cos 2 'y + ^ 0 /i, -f- J 7) 1; b 2 J t$bf»^ a » +|.J M-i- 1h ih ± j 

- . ‘ f ' , 

■ + ~ B x B a j cos 4 y + ( A t B, + *£A 2 ^ 2 ) sill 9 + (2/iji/ 2 -|- 

5- yLo //J sin 29 + (3^1 , B s + 2yl, B„ + yi 3 //,) slix 39. 


(T .26) 


( 7 . 27 ) 


( 7 - 28 ) 


( 7 . 29 ) 
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thje__^cinitX-.-Q.f . the rounded leading edge and the rounded trailing 
edge c~ <l 9 i.^ 1I jj , where relationship ( 7 . 25 ) already has subs- 

stantial de-viations., it .is appropriate to .calculate the contour velocity 
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using the formula: 






f ' ■ 


. 't 


&Z- 


_*• • . ; V ty, i* t ff/j 


Here, we have 


(7.30) 
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. fcjc — 1+ 5) (Z n feds7fi)- + fe— 3i+2(jf 0 + fflf)dori4>— . ttfecb&M&l- 4 ' *..) 

' - WZz ' 1 


■l.'t 


rt“U 

t t 4, 


!. -j- - - ■ '^j 0 s , m — j— •£ E, 's* Rin ^ + « -7— — --£ ^ 4 “ u j^? v ; (7.31) 

! : 1— feosS-4- TT^ 2 ,t '’ X ’ i" ~bOS^ L l-~fe 2 ' ' 5 1 


- * v/ „> »V \ 

2o , I u • 


' + 


■2o- 


— r —2{<f tf +ff 0 ).fcds , y' + fl h$ n QO&ii& 
*d--'’.' . *-2 j 


- l+cos&+^E- 

.6/*- '-,’V . 


s > 


- - ‘SoSbiS j -'1 {fosiii$ 1 . , , v, ’ . J 

X 7 f— ■ + — ■ £ j — 2.{ffo +So) Sin 9-f qn Hill n&— \ 

a — eos^-h— E s l-f bo&&+ 4-e? r • s. — • . »*“■ ••• 

< O * . ' O . 1 , * ' f 


-■:> q Q sin & . 


:*'*'•'* e 


C 'X 


t 2 
' it»V 


i ' — X fa- Cbs^)-l-'e f(fiJ 4- ^3i) + 25 o)' siii-9- + 

i‘ v ■■ »-t . h r . • - , r , . • ; - 

4 

+ ~ 4. 


j, 

, 1 


ffo ■ 


1 — bds-B-f — e* '"- 1 
2 


* / ' / 
2 n<j n cbs iiS-ji * 


The velocity along the leading edge is 


(7.32) 



(7.33) 


THE ODD PROBLEM. (THE FIRST MAIN PROBLEM) OF PROFILE THEORY 

The odd problem in the true sense of the word is to find the profile 
shape from a specified velocity distribution on the profile contour. 
Here., we will not 'solve this problem in this generality. We 'will in- 
troduce a modification so that the shape of the profile placed into 
the inhomogeneous primary flow field is calculated for a specified 
skeleton chord, source distribution and vortex distribution on the 
skeleton. 


The expression, for the components of the primary velocity at the 


40 


ORIGINAL page is 
OR POOR QUALITY 



point of the sought-after profile ils found by expanding the components- 
of the primary velocity in the vicinity of the skeleton chord in Taylor 
series. In this calculation, the jfirst- terms of this expansion are con- 
sidered only. Instead of (7.12), we will write 
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;.#r 

t 

f * - u 0 


y» + t cds‘ ria : 




(7-34) 


(7.35) 




,ji a '» { 

The coefficient values y«j v fc> .and ~T7: are 'determined from the 

. M. 

component values of the primary velocity and their derivatives at the 
chord points. 


The dash notation has been introduced in order to distinguish these 
coefficients y^ and from the coefficients y n and v n of the expansions 
of the primary velocity components; (7.21) calculated • directly on the 
profile. Using (7-3) and considering (7.2), we then have: 


& — 


% u itu 

! {*2 
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f 

: y-i 


V-i — 




S ^0 — 
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Expressions (7*36) and (7-37) are only used where the coefficients 
P n and v n appear in the first order terms. In the second order terms, 
we only write 

l^ri" ~ ‘f^H > ~ • (7.38) 


It should also be realized that the coefficients B n obtained by 
solving system (7.18), must satisfy condition (7.2). Therefore, every- 
where in system (7-18), we set: 
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(7.39) 


It follows from this that the coefficient j'gQ, when selecting the 
circulation distribution. Is not arbitrary, but instead is a result 
of solving system (7.18). 


After this preparation, we can begin with a solution of systems 
(7. 15)3 (7.18), (7.19). Prom the calculated coefficients A n , and 
C , we will then find the skeleton form using equation (7.1) and 
( 7 . 3)3 and after this the thickness distribution from equations (7*5) 
and ( 7 . 6 ). The total circulation and velocity distribution on the 
profile contour are calculated from equations (7.24) - ( 7 * 30 ), similar 
to the solution of the even problem. 


As already mentioned earlier, the velocity distribution on the pro- 
file contour is decisive for solving the odd problem. Therefore, we 
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will at least give an approximate expression for the contour velocity, 
corresponding to the selected vortex and source distribution. This is 
obtained from the expressions given above, by only considering the 
first-order terms. In a region outside the rounded trailing edges, /659 
the contour velocity is 
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The leading edge velocity is 
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The coefficient g Q is given by 


(/o' — v 0 


pi ! 

<J*l + Vgm ; 
4w 2 —1 . 


(7.41) 


(7.42) 


It is advantageous to estimate .the velocity distribution on the 
profile contour using the formulas (7.40)- - (-7.42)i, i n every case, 
before beginning the calculation, so as to exclude those circulation 
and source distributions which lead to profiles which have unfavorable 
aerodynamic characteristics. On the other hand, it is possible to use 
formulas (7-40) - (7.42) for selecting the coefficients g fi and q , at 
least in the first approximation when calculating the profile from the 
described velocity distribution on its contour. 


8. Examples 

The usefulness of the theory discussed in this paper will now be 
illustrated with three practical examples. 

8.1 Symmetric profile in a Homogeneous Plow Field 

As a first example, we will consider the computation for a modified 
NACA nrofile. The profile shape (Figure 11) is given by the coordinate 
table (Table 1) . 
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Figure 11 


Table 1 
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X = nondimensional leading -.edge distance 

i 

X = nondimensional profile ithickness 
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Using harmonic analysis, we determine the coefficients A n of the thick- 
ness distribution for this profile ,( see equation 7. 21 ) . 


4 - 


1 ~ 0,0753 ; , • ' Ai = 0^0192.; 




; ALj = 0;0002 \ • A b '= 6j0008 ; ' 


0,0065 ; 'j 
0,00(321 I 


! S nA * “ OjO&SzJ (-1)^ nl n ='~ 0,0226. 

J** ^”1 , «*"1 i * /. 

Since we have ^“syrmetric^prbTiXe^^ii 'o^f' the coefficients B (see 
equation (7.20) -are equal to 0. According to equation (7.19), the co- 

t 

efficients C n are: I 
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The velocities are referred to the x component of the homogeneous 

i ! 

primary velocity field so that according to equation (7-12), 
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where is the angle of attack. 


We will now introduce the angle 
rest of the calculation. All of t 
.pansions (7.12) are zero, with the 


of attack parameter K = tg a ffi for the 
he coefficients p and v in the ex- 
exception of v Q , for which we have': 


{ 

i 


These values are then substituted in equations (7 • 15)-(’7 . 1-8) , and /66l 
after this they are solved for the 'desired vortex and source distri- 

P 

bution coefficients I 


. % = 1,093 if • g'l=— 
, jo = OjC)5M‘;.g6 = — 


6 , 0 ' 40 ! k ] r g 2 = 0,oIik± Us '=■ 
Qfoiol ; q 2 = 'dj 0344 q' & = ij'C 






According to (7*240, the total circulation is 

Vicki = 3,371 7 -K. j 

i 

-The ve-locity distribution on the profile contour is obtained by sub- . 
stitution in equations (7.25) - (7- I 33) J and is shown in Figure 12 for 
the angles of attack a ro = 0 and two degrees. The velocity distribution 


I 
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8.2 Curved Profile in a Homogeneous Flow Field 


As a second example, we will now consider .a profile with curvature 

I 

in a flow, which was produced by plotting a symmetric profile mentioned 
above on a circular arc profile center line. The central angle of the' 
arc shaped center line is w = 2.9°32’. 

f /662 

i . . 

For circular arc center lines, the coefficients in expansion (7.20) 

are given by the following expressions: 
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In the example given, only two coefficients B n are different from zero: 


B 1 = 0.2613 


B^ = 0.0022 
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The total circulation is 

7 r 


, ~ (J;4690 + 3,371 K ! 

■JyU'* ; . -j 


/6 


and the velocity distribution on the profile contour is shown in Figure 
13 , for = 0°26', which corresponds to an entry without shocks. We 
also show it for a = - 1°34’ and 2°26'. 


i 

8.3 Curve Profile in Cascade Configuration 

i 

Finally, we will calculate the flow around a profile from a previous 
example, with a cascade configuration (Figure l4). The division ratio 



Figure l4 

i 

i 

(s/c) and the separation angle (A) of this cascade are given by 

s/c = 0.904 A = 34°30 ! 

The flow around the profile in a cascade configuration is considered 
as the flow around an isolated profile in a nonhomogeneous primary flow 
field [7]* The inhomogeneity of the primary flow (induced velocities) 
is caused by the other ’ cascade profiles. 

i 

... ... I. 

The x component (U .) of the -translation- velocity (W ) is selected 

I 

a-s the reference velocity. The translation velocity is defined as the 
vector average of the incident and; departing flow velocities (W^, W 2 ). 
The coefficients and v n in the expansions of equation (7.12) are 
given by -the following expressions: (see [8]): 
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Figure 15 
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The calculation is completely similar after this. Of the numerical 
results , we will only give the departing flow angle a p (Figure 15) 
and the velocity distribution on the profile contour. The velocity 
distribution is shown in Figure l(j for entry without shocks (a^ = ll 0 i33') 
and also for ct^ = 9°33 T and 13 0 33'|* Figure 16 shows the velocity dis- 
tribution referred to the incident- flow velocity. 


CONCLUSION j 

This paper contains a theory of a moderately curved profile in a 
nonhomogeneous f 1 o.w_f.ie.l JTh e~t n e or y - Asihuiit.~up.~as a "second order 
theory". The coefficients in the expansions for the circulation source 
and thickness distributions, as wejll as the coefficients in the expan- 
sion of the skeleton shape, are considered !asj small "of first order". 
All third-order quantities are ignored in the relationships . 

The main results of the paper are adapted to the requirements for 
numerical calculations in the seventh chapter.' Without knowing the 

.theory in detail, it can be used J This circumstance is very impojr- 

tant , because the derivation of the theory requires a knowledge of 
functions of a complex variable, and conformal mapping. In a numerical 
calculation, all that is required |is the knowledge of expansions of 
trigonometric series (harmonic analysis) and the solution of systems 
of linear ' inhomogeneous equations. This assumption can be assumed 
to be satisfied for any technical employee. The systems of linear 
equations have a very convenient fjorm for numerical calculation. The 
largest quantities are always located along the main diagonal. This 


means that we can recommend a step-wise approximation method for the j 
solution, and in this case it also leads to a rapid solution, even for 
several unknowns. } 


In practice, we often encounter this type of flow, and this theory 
-was^ -formulated in general terms without direct adaptation to spscAal^ 5 
cases. We can consider -the following applications: flow around aero- 4 
dynamic profiles in a homogeneous jflow, flow around profiles near the 3 
earth, flow around straight blade cascades, flow around radial blade 2 
cascades, flow around guide bladesj, etc. In all these cases, we can 1 
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transform this theory into calculation formulas which will consider the 
specific properties of the individual applications. This results in a 
further simplification of the calculations. Of the technical applica- 
tions, the calculation of the flow around straight and radial blade 
cascades is probably the most important. Consequently, the results of /666 
this paper and the results of papers [ 7 ] and [8] were formulated into a 
calculation procedure for calculating the flow around blade cascades 
with small thicknesses and moderately curved blades. This has been done 
for incompressible and compressible subsonic flow, as already discussed 
in example ( 8 . 3 ) . 

The formulas for the numerical calculation discussed in chapter 7 
can be performed with an .electric calculator or a slide rule. They 
are very well suited for programming of automatic computer installations. 
The time requirement for a complete calculation of a profile or a blade 
cascade is reduced to a minimum.. In Chapter 8, we discussed examples 
which have already been calculated on the ZUSE-23 computer. 

The results can be generalized to profiles with a high curvature, 
but with a rather small thickness, as has been shown in paper [6] for 
the flow around a thin, highly-curved profile in a nonhomo geneous flow 
field. These generalizations were not included in this paper, because 
moderately curved profiles are.'very important in technical applications. 

Also inclusion of this general case would have made this paper much more 
cumbersome. 
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